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Abstract. A common challenge to proving asymptotic stability of solitary 
waves is understanding the spectrum of the operator associated with the lin- 
earized flow. The existence of eigenvalues can inhibit the dispersive estimates 
key to proving stability. Following the work of Marzuola &; Simpson, we prove 
the absence of embedded eigenvalues for a collection of nonlinear Schrodinger 
equations, including some one and three dimensional supercritical equations, 
and the three dimensional cubic-quintic equation. Our results also rule out 
nonzero eigenvalues within the spectral gap and, in 3D, endpoint resonances. 

The proof is computer assisted as it depends on the sign of certain inner 
products which do not readily admit analytic representations. Our source code 
is available for verification at htt p : //www. math. toronto .edu/slmpson/flles/l 
|spec_prop_asad_simpson_code .zip! 



1. Introduction 
The nonhnear Schrodinger equation (NLS) in M'^+^ dimensions, 

(1.1) z7/.t + VV + /(l^nV' = 0, V(0,x) = ^o(x), 

appears as a leading order approximation to a wide variety of natural and engineered 
systems, including optics, plasma physics, and fluid mechanics, ft is also a model 
equation for the studying the competition between nonlinear and dispersive effects. 



See Sulem & Sulem, 21 , for details and examples. 



1.1. Solitons. For appropriate choices of the nonlinearity, / : M — > M, the equation 



will possess solitary wave solutions, smooth localized solutions of (1.1 1 satisfying 
the ansatz 

(1.2) x) = e*"*i?(x; w), uj > 0. 
R then solves the nonlinear elliptic equation 

(1.3) -ujR + V^R + f{\R\^)R^O. 

The solutions are assumed to have finite norm and vanish at infinity. 

These nonlinear bound states, deemed solitary waves or solitons, are of interest 
both as mathematical objects and in applications. This raises the question of 



whether or not these solutions are stable to the flow of (1.1). Developing tools for 
assessing stability is the main purpose of this paper. 
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1.2. Stability. With regard to solitary waves, the two types of stability typically 
discussed are orbital stability and asymptotic stability. Orbital stability, alterna- 
tively called Lyapunov or modulational stability, asserts that if the data is close to 
a solitary wave, it remains close, modulo a group of symmetries associated with the 
equation. 

In [22"23', M.I. Weinstein proved that orbital stability of a solution to ( 1.3 1 could 
be assessed by computing the sign of 

(1.4) A y'|i?(x;c.)|2dx. 

When it is positive, the solitons are stable; when it is negative, they are unstable. 
This work was subsequently generalized for Hamiltonian equations by Grillakis, 
Shatah, & Strauss [9)[To] 



For NLS with a power, or monomial, nonlinearity of the form /(s) = s'^, the 
transition between stable and unstable solitons is determined by the product ad. If 
ad < 2, the solitary waves are stable, otherwise they are unstable. These regimes, 
crd < 2, = 2, > 2 correspond to the subcritical, critical, and supercritical forms of 



(1.1), and are intimately related to the well-posedness of NLS, |21) . 

Of course, orbital stability does not tell us of the limiting behavior of the solution. 
An orbitally stable soliton could perpetually oscillate through the symmetry group 
of the equation. It is our expectation that perturbations of stable solitons diminish 
as t — >■ oo, and the solution relaxes to a particular solitary wave with a fixed set of 
parameters. We turn to asymptotic stability to understand the dynamics as i — >■ oo. 



Asymptotic stability is usually proven by expanding ( 1.1 ) about the solitary wave 
solution, ip — e"^*(i?(-; w) + u + iv) to arrive at an equation for the perturbation, 

p — u + iv, 

(1-5) 

^ 1\ fL+ \ fu 



1 OJ \0 L^J \v 

F contains terms that are nonlinear in the perturbation, and the scalar operators 
L± are given by 

(1.6a) L+ = -V^ +LU- f{R^) - 2/'(i?2)i?2 ^ _y2 + ^ + 

(1.6b) = -V^ +UJ - f{R^) = -V^ +LO + V-. 

Proving asymptotic stability of the solitary wave will require Strichartz estimates 
of the form 

(1-7) lle^^VLrLj <ll/IUj. 

With such an estimate, one can conclude linear stability of the soliton from the 
temporal decay of a solution to the linear problem 

Pt = JLp, p = (u, w)^. 

This decay can then be used to show that the nonlinear part of the flow, F in 



(1.5), is dominated by the linear part can can be treated perturbatively. Successful 
implementations include Buslaev & Perelman, [T|, Buslaev & Sulem, [2|, Cuccagna 
and Rodnianski, Soffer & Schlag, [SlfTs] and more recently to Schlag and Krieger 



& Schlag, 11 16 . These last two works, where the authors show the asymptotic 
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Figure 1. The generic spectrum of JL for an orbitally unstable 
soliton. The instabihty is due to the positive eigenvalue, /j,^. 



stability of a constrained soliton which is orbitally unstable, are closely related to 
the present results. 



1.3. The Spectrum & Embedded Eigenvalues. Estimates of the form (1.7) 
require adequate knowledge of the spectrum of JL, a{JL). Since the solitons are 
highly localized, JL is easily shown to be a relatively compact perturbation of 

-V^ + 

y^-uj 



(1.8) 



which has as its essential spectrum 

{—ioo, ~iio] U [iu:,ioo) 
This can easily be computed by the Fourier transform. Thus 
(1.9) (JcssiJL) ~ i(— oo, — w] U i[u!, oo) C a{JL) 

See [t] for additional details. By direct computations, one can show that the origin 



is an eigenvalue of JL of algebraic multiplicity at least 2d + 2. See section 2.1 below 
for the elements of the kernel, and Figure [l] for a visualization of the spectrum of 
the problems considered here. 

The existence of purely imaginary eigenvalues both within the spectral gap, 
[—iujyiuj], and embedded in the essential spectrum and resonances can obstruct 
the Strichartz estimates. Many results on asymptotic stability of NLS solitons. 



including [2||7||16j, have explicitly assumed: 

(1) JL has no embedded eigenvalues; 

(2) The only eigenvalue within the spectral gap is zero; 

(3) The endpoints of the essential spectrum, zticj, are not resonances. 

We are thus motivated to find ways of rigorously proving these assumptions. 
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Some previous results on these spectral questions were addressed by Krieger 
& Schlag in ID for /(s) — s" with > 4 (supercritical), [TT]. Alternatively, 
Cuccagna & Pelinovksy and Cuccagna, Pelinovsky, & Vougalter [4j[5] showed that 
in ID, embedded eigenvalues with positive Krein signature can exist and not disrupt 
asymptotic stability. But an extension of their approach, based on the Fermi Golden 
Rule, to higher dimensions remains a challenge; it may be easier to prove that such 
states are absent. In [6], Demanet & Schlag gave a numerically assisted proof for 
the absence of non zero eigenvalues within the spectral gap for 3D NLS with a 
power nonlinearity with 

(1.10) 0.913958905 ±l(f <a<l 

More recently, Marzuola & Simpson, |l2j, proved that the 3D cubic problem 
has no purely imaginary eigenvalues or endpoint resonances. This proof required a 
modest amount of numerical assistance to compute: 

• The dimension of the negative subspace of several ID Schrodinger opera- 
tors; 

• The inner products of solutions of several ID boundary value problems. 



In this work, we apply the approach of 12 , to some other equations 



1.4. Main Results. Our main results are for the ID and 3D supercritical NLS 
equations with /(s) = s" ^ and the 3D cubic-quintic nonlinear Schrodinger equation 
(CQNLS) and for the ID supercritical NLS, 

(1.11) zV't + VV + (lV'l'-7l^l^)V^ = 0. 
Theorem 1. For 3D NLS with 

(1.12) 0.807425 < a < 1.12092 

JL, the linearization about the soliton R — i?(-; 1), has no purely imaginary eigen- 
values or endpoint resonances. 

Theorem 2. For 3D CQNLS with 

(1.13) 7 < 0.00989115 

JL, the linearization about the soliton R = R{-; 1), has no purely imaginary eigen- 
values or endpoint resonances. 

For the sake of exploring the limits of our approach, we also prove: 

Theorem 3. For ID NLS with 

(1.14) 2.4537956056 <a < 6.1288520139 

JL, the linearization about the soliton R — R{-; 1), has no purely imaginary eigen- 
values. 

As previously mentioned. Theorem [3] has been rigorously established for all a > 
2, ~ 
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To prove these theorems, we first show that a particular bilinear form, which 
would vanish at any purely imaginary eigenvalue or endpoint resonance, is coercive 
on an appropriate subspace. This coercivity allows us to immediately rule out such 
states. To proceed we must define distorted variants of L± which form the bilinear 
form. 
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Definition 1.1. Given L± and a skew adjoint operator A, define the two Schrddinger 
operators via the commutator relations: 

(1.15) C±^^[L±,A]=-V' +V±. 

For z = {u, v)^ E X , define the bilinear form 

Biz,z)=B+{u,u)+B-iv,v) 
— [L+u, u) + {L-v, V) 

The operator JL is said to satisfy the spectral property on the subspace U C x 
if 

(1.17) S(z,z)> J (|Vz|2 + e-l"l|z|2)dx 

The skew adjoint operator A and the subspace U are, at this point, unspecified. 
In our work, we use 

(1.18) A=^+x-V. 
The potentials, V±, this induces from L± are 

(1.19a) V+ = ix • V [/(i?2) + 2/'(i?2)i?2] = r(3/'(i?2) + 2f" {R^)R^)RR' 

(1.19b) V_ = ix • V [/(i?^)] = rf{R^)RR 

For the 3D problems we consider, these correspond to: 

f{s) = s" :V+ = (7(2(7 + l)i?2'""^i?', V_ = arR^^-^R' 
f{s) = s- :V+ = r(3i? - 10-/R^)R', V_ = r{R - 2-/R^)R' 



The particulars of this subspace will be discussed in Section |2.2| This Spectral 
Property appeared in the works of Merle & Raphael, and Fibich, Merle k Raphael 
in their proofs of the log — log blowup of critical NLS, [S 13 . It has similarly 



appeared in Simpson & Zwiers, in a proof of the log — log blow up for vortex solitons 



in 2D cubic NLS, 20 



The role of this Spectral Property was identified by G. Perelman, [14|, and leads 
directly to: 

Theorem 4. // the Spectral Property holds for JL, then JL has no imaginary 
eigenvalues on hi. 



The proof is quite simple and appears in 12 . Briefly, one shows by direct com- 
putation that B vanishes at any eigenstates corresponding to imaginary eigenvalues. 
In 3D, this same analysis will rule out endpoint resonances. 

Though all of the cases for which we have results correspond to orbitally un- 
stable solitons, the results remain of interest. They can be used to prove results 
on constrained asymptotic stability, as in 16 , where one projects away from the 
linearly unstable direction of JL. In addition they are intrinsically interesting as 
results on the spectrum of Hamiltonian operators. Finally, these theorems also map 
out the scope of success for this approach. 

Our paper is organized as follows. In Section[2]we review some relevant algebraic 
properties of JL, and define the subspace, U, that will be used to prove the Spectral 
Property. In Section [3j we present some preliminary results used in the proof. 
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Sectiorj4] presents our computations and the main resuhs. Some discussion of the 
hmits of this approach are raised in Section [5] 

Acknowledgements: The authors wish to thank J.L. Marzuola for some helpful 
comments, D.E. Pelinovsky for suggesting a comparison with the Demanet & Schlag 
threshold, and C. Sulem for suggesting the extension to CQNLS. This work was 
supported in part by NSERC. 

2. Algebraic Structure 

As noted in the introduction, a{JL) includes the essential spectrum, lying on 
a portion of the imaginary axis, along with some eigenvalues. Since the algebraic 
structure of JL is used in our proof, we review the discrete spectrum here. 



2.1. Discrete Spectrum. Given that (1.1) supports a solitary wave solution, we 
readily observe: 

Theorem 5 (Kernel). JL has a kernel of algebraic multiplicity at least 2d + 2: 

forj = l,...d. 

For power nonlinearities, we can, and do, use, that 

(2.1) 2d^\^^iR= ^R + x-VR. 

In addition to these elements, there are two additional eigenstates that can reside 
on the imaginary axis, the origin, or the real axis, depending on the slope condition. 



(1.4). For orbitally unstable problems we have: 

Theorem 6 (Off Axis Eigenvalues). There exists a pair of real, nonzero, eigenval- 
ues located at ±/i*. Corresponding to /i^, > 0, is the eigenstate cf) = (0i,02)"^- For 
d = 1, the 4>j are even functions, and for d > I, the (f>j are radially symmetric. 

These off axis eigenvalues appear in Figure [T] Their existence follows form the 
results of Alternatively, Schlag presented a continuity argument in ^16J that 
relies on the knowledge of the two additional eigenstates that appear at the origin 
in the critical case, ad — 2. 

It is well known that for f{s) = s"' , the slope is negative when ad > 2. Due to 
scaling, this holds for all u > 0. For CQNLS, some solitons are stable and others 
are unstable, depending on 7 and the soliton parameter uj. In this work, we shall 
fix w = 1 for which solitons arc known to exist provided 7 < 3/16, |21 . We now 
determine the range of < 7 < 3/16 for which the solitons are unstablej^ 

At 7 = 0, this corresponds to the 3D cubic problem, which we know to be un- 
stable. By continuity, we expect there is some open set of 7 near zero for which the 
solitons are orbitally unstable. To identify the threshold value of 7, we numerically 



compute (1.4 1, and find that it changes sign at 
(2.2) 7* « 0.0255453. 

The slope condition values as a function of 7 are plotted in Figure [2] See Appendix 
[K\ for details of this computation. 

^Alternatively, we could set 7=1, and vary oj. 
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Figure 2. The numerically computed ( 1.4 1 at w = 1 changes from 
negative to positive at the critical 7+ ~ 0.0255453. 



2.2. Subspaces. As discussed in Section 2.2 of [12], th e subspace U for which we 
wish to establish the Spectral Property of Definition |l . l| must satisfy two properties: 

• If they exist, eigenstates of purely imaginary eigenvalues must reside in U] 

• U must be orthogonal to the negative subspaces of C±. 

Obviously, if both properties are satisfied, then there are no purely imaginary eigen- 
values. But it may be possible to construct a U that is orthogonal to the eigenstates 
of such eigenvalues, and on which the operators are positive. 

A practical subspace makes use of the eigenstates of the adjoint operator, (JL)* . 
These are closely related to the elements of JL. Given i/; = (-01, -02)^, an eigenstate 
of JL with eigenvalue A, 

(2.3) iJLr(tA^-xf^^ 



ipij yip 



Consequently, 



Proposition 2.1. Let xjj ~ (01,02)"^ be an eigenstate corresponding to a purely 
imaginary eigenvalue of JL. 

• Any element ofkeTgdJL)*) is orthogonal to tp; 

• The unstable eigenstate (f> of JL satisfies the relations: 

(Vl,02)=O, (V2,0l)=O. 



See Section 2.2 of 12 for a proof. Using this, we define U to be the orthogonal 
complement of 

-{©•u".)'(.:«)'(T).eo).(;. 

Indeed, 

Theorem 7. If the Spectral Property can be established onlA given by the orthogonal 



complement to (2.4), or any subspace containing lA , JL will not have any purely 
imaginary eigenvalues on L^ x L^ . 

We succeed in proving: 
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Theorem 8. Assume z = {,f,g)'^ £ x satisfies the orthogonality conditions 
(2.5) {f,R) = (/,02) = (.9,a^U=ii?> = (5,0i> = (/,xi?) =0 

Then the Spectral Property holds for: 
3D NLS: Provided 

0.807425 <a < 1.12092; 

3D CQNLS: Provided 

7 < 0.00989115; 

ID NLS: Provided 

2.4537956056 < a < 6.1288520139. 



(3.1) 



(3.2) 



For ID and 3D NLS with power nonlinearities, we use (2.1). 

3. Preliminaries 

We shall prove the coercivity of B in three steps, closely following f8','T2', and 
omitting some details. 

• In ID, we decompose into even and odd functions and study 

B+ {fj)+ B. {g,g)= B^ {&\f'^''^)+ B^°^ (/(°) , /(°) ) 
+ B^l\g^^^\g'^^y)+B^^\g'^°\g^°y) 

We then show that on U, each of the four forms, B^^^°\ is positive. 
Similarly, in 3D, we decompose into spherical harmonics, 

oo Lfc 

B+ if, f) + 6_ (5, ,g) = E E (/^"'^ ' /'""^ ) 

CO Lfc 
k=0 1=0 

where /C^'') and (7^^'') are the components of / and g in harmonic k with 
angular component I. As C± are radially symmetric operators, they are 

(k) 

independent of I. Here too, we shall prove positivity of each B^. onU. 

• For each of these bilinear forms, we shall determine the codimension of 
a subspace on which they are positive, called its index. This is problem 
specific and can vary with the dimension and choice of nonlinearity. 

• Lastly, we shall show that the orthogonality of / and g to U-^ induces 
orthogonality, with respect to to the bilinear form, to the negative subspaces 
of the form. See Section 4.1.3 for an example. 

3.1. Indexes and Eigenvalues. The index of a bilinear form S on a vector space 
V is defined as 

indy(i3) = max{/c e N jthere exists a subspace V C V 

of codimension k such that -B|y is positive} 

To compute the indexes, we rely on the following propositions: 
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Proposition 3.1. Let U^''^ and U^"'' solve the initial value problems 

LC/(«) =-^u(«) + y(a:)i7(^) =0, U^^\0) = 1, ^f/(«)(0) = 

LC/W =--f!^u('') + l/(x)C7(°) =0, ?7(°)(0) = 0, -^U^"\0) = 1 
ax'' ax 

where V is even, sufficiently smooth, and \V{x)\ < e"**'^', k > 0, and let 

Then the number of zeros of U^'^^ and U^°^ is finite, and 

indiifi(-B) = number of positive roots ofU^'^^ 

= number of negative eigenvalues in 

indjji (B) = numher of positive roots of U^°^ 

= number of negative eigenvalues in 

Hi and are the even and odd subspaces of H^. 

Proposition 3.2. Let U^^^ solve the initial value problem 



l{k)u(k) ^ 



dr^ r dr 



C/(^) = 



lim T-^U^^\r) = 1, lim — ( r'^U^'^Hr)] = 0, 

r-)-0 ^ ' r-)-0 dr \ J 

for A: = 0, 1, 2, 3, . . . where V is sufficiently smooth, and \V{r)\ < e~'^l'"l,«; > 0, and 
let 

Then the number of zeros of U^^^ is finite, and 

indjji ^(-B^°^) = number of positive roots ofU^"'' 

= number of negative eigenvalues in HI^^ 
indjji ^^{B^^^) = number of positive roots of U^''\ for k = 1,2, . . . 
= number of negative eigenvalues in H^g^^ 
The space -ffra,d+ subspace of radially symmetric functions for which 



/ 



i/hxr'rfx<oo. 



The indexes in the different harmonics have an invaluable monotonicity property 
that permits us to restrict our attention to a finite number of harmonics: 

Corollary 3.3. Fixing the potential V{r), the bilinear forms B^''^ satisfy 

ind(B(*^+i)) < ind(B('=)) 

The last result that we shall make use of in computing the index of a bilinear 
form is that it is stable to perturbation by a sufficiently localized potential: 
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Proposition 3.4. Fixing an operator L from Propositions\3 . j] or\3//i\ or Corollary 



3.3. there exists 5q sufficiently small such that the bilinear form, B, induced by the 
perturbed operator 

satisfies 

ind(B) = ind(B) 

The proofs of the preceding results are given in fl^, and references therein. As 
one might suspect, this is closely related to Sturm oscillation theory. 

3.2. Boundary Value Problems. In the course of proving the spectral property, 
we will need to solve a series of problems of the form 

Cu = f 

for u where C is one of the above operators and / is a localized. Though the 
solutions of these problems are found numerically, the invertibility of the operators 
can be rigorously justified, up to the index computation. For the index, we continue 
to rely on numerics. 

Proposition 3.5. Let f be a smooth function, with even or odd symmetry, satis- 
fying the bound |/(a;)| < e~'''^' for some k > 0. Then there exists a unique solution 

ueL°°(M)nc2(M) 

to 

Cu = f. 

where C is one of C±. u possesses the same symmetry as f. 

Proposition 3.6. Let f be a smooth, radially symmetric, function satisfying the 
bound \f{r)\ < e"**^ for some k> 0. Then there exists a unique solution 

(1 + r''+^)u e L°°{R^) n C^{R^) 

to 

Cu^f 

where L is one of C^^^ . u is radially symmetric. 

While the solutions of the 3D problems decay cx r^^^*^, the solutions of the ID 
problems are merely bounded. Knowledge of the index of each of the operators is 
necessary in proving the uniqueness of the solutions. Additionally the invertibility 
of the operators is stable to perturbation: 



Proposition 3.7. For sufficiently small 5q > 0, the results of Propositions 3.5 and 
\3.()\ will also apply to 

£ = £ - (Joe-I'^l 
4. Results 

We are now in the position to present our results and prove Theorem [8] for each 
NLS problem. From this, we can conclude Theorems [l] [2j and |3] Throughout, we 
shall make use of the orthogonality conditions (2.5 1. See Appendix [A| for details of 
our numerical methods. 



4.1. 3D Supercritical NLS. 
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(a) Harmonic A; = (b) Harmonic fc = 1 




1.4I ' ' ' ' ' 
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(c) Harmonic k = 2 

Figure 3. Index functions for the 3D supercritical NLS equation 
computed at six values of .75 < a < 1.5, inclusive. 



ind(e^"') = ind(Si"^) = ind(S^'^) = 1 



4.1.1. Indexes. 

Proposition 4.1. For .75 < a < 1.5, 

ind(6^^^) =ind(Si^^) = 

Proof. Examining Figures [3j we see that for six computed values of .75 < ct < 1.5, 
inclusive, we have the corresponding number of zero crossings and apply Proposition 
3.2[ We argue by continuity that this should hold at all points in the interval. These 
were computed using the approach described in Appendix ] A. 1[ □ 



Corollary 3.3 ensures positivity of the higher harmonic forms. 



4.1.2. Inner Products. 

Proposition 4.2. Let Jj}"-* and U2^\ solve 

(4.1a) C^^^U["^=R 

(4.1b) /:fc/f=02 
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And define: 
(4.2a) 
(4.2b) 
(4.2c) 

The Kf^ have the values indicated in Figure^^ Moreover, for .8 < a < 1.2, 

(4.3) iK['^Ki'^ - {Ki'^f)/4"' < 

as pictured in Figure 



iff 

iff L/(°\c/f ). 



Proof. Proposition |3.6| ensures that the solutions exist. The boundary value prob- 
lems and inner products were computed using the techniques described in Appendix 
|A) The (|42| and (|43| were then computed at 41 uniformly spaced values of a be- 
tween .8 and 1.2. □ 
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Figure 5. (4.3 1 as a function of a for 3D supercritical NLS. 



Proposition 4.3. Let zf'^ and Z. 



7(0) 



live 



(0) 7(0) 



t(0) _ 
Jl — 



(4.4a) 
(4.4b) 
and define: 
(4.5a) 

(4.5b) 

(4.5c) 

The jj*'-' have the values indicated in Figure Moreover there exist ci > .8 and 
(72 > -8 such that 



J. 



(0) 



£(0)^(0) ^(0) 

/-(0)7(0) 7(0) 
L,_ , Z.2 



(0)^(0) 7_(o) 



jf ^ < 0, /or CTi < cr < 1.2 
/f)^)/j( 



(4.6) 
(4.7) 

as pictured in Figure where 

(4.8a) cTi = 0.807699 

(4.8b) (72 = 0.807425 

Proof. Like in the proof of Proposition |4.2[ we numerically solve the boundary value 
problems at 41 uniformly spaced values of a between .8 and 1.2. Using the root 
finding approach, discussed in Appendix |A.3[ we found the zero crossings, ai and 
(72, of j'^i'^ and (??) to be at the indicated values. □ 



Proposition 4.4. Let u'l^^ solve 

(4.9) 
Define 
(4.10) 



rR, C/}^^ e L° 




Figure 7. (4.7| as a function of a for 3D supercritical NLS. 



EIGENVALUES AND NLS 



15 




Figure 8. k[^'' as a function of a for 3D supercritical NLS. 



There exists > .8 such that 

k[^^ < 0, for .8 < CT < CT3 
as pictured in Figure where 
(4.11) 0-3 = 1.12092 



Proof. This is solved in the same manner as Propositions 4.2 and 4.3 



□ 



4.1.3. Proof of the Spectral Property. Subject to the acceptance of these computa- 
tions, we are now in the position to prove the spectral property for 

0-2 < CT < CTa 



where CT2 and era are as defined by (4.8a) and (4.11). We first shall show that the 



bilinear form on each harmonic, b''^\ is positive. We first prove the case of b'^^K 



Fixing cr, for a sufficiently small 6o, in the sense of Propositions 3.4 and 3.7 
1 



K. 



(0) 



>(0)n2 



< 0. 



where A"j°^ is the inner product associated with the perturbed operator, 

Given that / is orthogonal to R and 02, it is orthogonal to any linear combina- 
tion, including 



R 3 



(0) 



K. 



(0) 



Let Q solve 
(4.12) 

By construction, 

(4.13) Sf(Q,Q) 



(0)^(0) 



(0)^2 



< 0. 



To complete the proof of positivity, suppose Q G -ffjadi even thought it is not; it 
decays too slowly for L"^. We could decompose H^^^^ as 

-f^rad = span {Q} © span {Q}^ 
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where the orthogonal decomposition is done with respect to . We can do this 

because ( |4.13[ ) imphes the form is non-degenerate. Since ind(S'^'') = 1, > 

on span {Q} . Were this not the case, it would yield a second negative direction, 
independent of Q, contradicting the index computation. 

Given / £ i?4d' / -L ^ and / _L 02j with respect to orthogonality, / is also 
orthogonal to q. We might then decompose / as proposed in the preceding 
paragraph 

f = cQ + f^, ^f(/^,0) = 0. 

If c = 0, then f = f'^ resides in a subspace of -ff^ad where B^"^ > , completing 
the proof. We now show c = 0. 

Taking the inner product of / and q, 

= c (g, g) + {f^,q) - cSf (Q, Q) + (^^, £f q) 

= cBf (Q, Q) + Bf {f^,Q) = cBf (Q, Q) + 0. 

Since we have computed B^^\Q,Q) ^ 0, c = 0. Hence, orthogonality to q 
induces orthogonality to Q. Thus, / hes in the positive subspace of 

However, Q is not in L^. To make this argument work, one can introduce a 
cutoff function and take an appropriate limit. We omit these details and refer the 
reader to [8 12 



Positivity of B^^ is proven similarly. In this case, we use the orthogonality 
of (7 to ■^R + tR' and (t>i. Indeed, we construct a new Q, the solution to 

■h 

However, this is only successful for (T2 < c < 1.2. 

Positivity of b'^^^ is somewhat easier, as there is no need to form a linear com- 
bination of elements; there is only one direction, ri?, to project away from in this 
harmonic, k'^^ < for all the values of .8 < cr < CT3. We conclude that orthog- 
onality of / to rR yields positivity. Since all other forms have index zero, there is 
nothing to prove for them. 

We have now established that 

(4.14) g(z,z)>0 

since each form on each harmonic is positive. It is trivial to see that 

(4.15) B{z,z) > So J e-l^l|z|2dx. 
This is almost the desired expression. Consider, 

(4.16) 6(z,z) > (^j |Vzp - IV+II/P - |V_||.gpdx^ + (1 - 6)5o J e-^-^\z\^d^ 

for all 9 G (0,1). Both potentials satisfy the estimate |V±| < e"^'^'''', cr > ^, as 
|x| — > 00, so taking sufficiently small, 

(4.17) 6(z, ^)>oJ iVzpdx + J e-|x| j^p^^. 
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Hence, 



B(z,z)> J (|Vz|2 + e-l"l|z|2)dx 



This proves the Spectral Property from which we then immediately get Theorem 
[T] for purely imaginary eigenvalues. As resonances in d = 3 have sufficient decay, 
we can also rule them out. 

4.2. 3D CQNLS. In this section we prove the Spectral Property for the 3D 
CQNLS equation. It is quite similar to 3D NLS, though we are now concerned 



with the values of 7 in (1.11 1 for which it holds. 

4.2.1. Indexes. 

Proposition 4.5. For < 7 < .012, 

■mdCf ind£(°^ = ind£i = 1 

indC^^'' = mdC^_l^ = 

Proof. Examining Figures |9j we see that for three computed values of < 7 < 
.012, inclusive, we have the corresponding number of zero crossings. We argue by 
continuity that this should hold at all points in the interval. □ 

4.2.2. Inner Products. 

Proposition 4.6. Let J/}"-* and U2^\ solve 
(4.18a) C^_^^U["^=R 
(4.18b) Cfu^°^ = 02 

Define xf^ as in ( |4.2[ ) . The k''^^ have the values indicated in Figure 
forO<-f < 0.012, 



10 Moreover, 



(4.19) (if("'if(°) - (x("))2)/if(°) < 

pictured in Figure [77} 



as 

Proof. As before, we prove this by direct computation, at twenty five values of 7 
between and .012, inclusive. □ 

Proposition 4.7. Let z[^^ and z'f\ solve 
(4.20a) = 
(4.20b) ^ - 



Define jf as in (4.5 I. The j'f^ have the values indicated in Figure 



there exist 71 > and (72 > such that 

(4.21) ^ < 0, /or > 7 < 71 

(4.22) (Jf - (Jf )^)/jf < 0, /orO<7<72 
as pictured in Figure where 

(4.23a) 71 = 0.00989115 

(4.23b) 72 = 0.0109065 
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1 2 3 4 5 

r 




■ 1 2 3 4 5 

(c) Harmonic k = 2 



Figure 9. Index functions for the 3D CQNLS equation computed 
at three values of < 7 < .012, inclusive. 



Proposition 4.8. Let U^^ solve 

(4.24) 



£WC/I^)=ri?, [/Wei- 



Define k[^' as in (4.101. Then 



iff ^ < 0, for < 7O.I2 



as pictured in Figure I4 



4.2.3. Proof of the Spectral Property. The proof is quite similar to that of 3D NLS 
and we omit many of the details. For a fixed 7 within the allowable range, we take 
5o sufficiently small so as not to alter the indexes or appreciably change the inner 
products. 

A notable difference is that for we let Q solve 

^(0) 

(4.25) 40)Q = g-=__|_i? + <^2. 
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0.002 0.004 0.006 0.008 0.01 0.012 

7 



(c) 

Figure 12. The jj"^ as functions of 7 for 3D CQNLS. 




0.005 0.01 

7 



Figure 13. {j[ '4 - (A ?)/J2 as a function of 7 for 3D CQNLS. 
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0.002 0.004 0.006 0.008 0.01 0.012 

7 



Figure 14. /sT^ ^ as a function of 7 for 3D CQNLS. 



By construction and Proposition 4.6 



(4.26) 



(Q,0) = 



1 



K 



(0) 



K. 



(0) 



< 0. 



If / Is orthogonal to R and 02, then it is orthogonal to any linear combination. 
As before, this will induce ,8^ orthogonality to Q, ensuring that S^'' is positive for 
such /. This holds for all < 7 < 0.012. 

Positivity of B^^^ is proven as before, except it only holds for < 7 < 72 due 



to the results of Proposition 



4.7 



Positivity of B^}'^ is the same, as are the rest of 



the forms, since they have index zero. This proves positivity of B on lA, and the 
coercivity oiB oxiU in the form of (1.17) follows as before. This yields Theorem [2] 



4.3. ID Supercritical NLS. In contrast to the 3D problems, where we decom- 
posed into spherical harmonics, in ID we decompose into even and odd functions. 

4.3.1. Indexes. 

Proposition 4.9. For 2.3 < ct < 6.3, 



ind(S^'^) = ind(si"') = ind(6^°^) = 1, 



ind(6i°') = 



Proof. As before, we prove this by direct computation. The index functions appear 
in Figure [T5| □ 



4.3.2. Inner Products. 

Proposition 4.10. Let uf''^ and solve 



(4.27a) 
(4.27b) 
and define: 
(4.28a) 
(4.28b) 



£(e)^(e) 



R 



C 



K. 
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X X 



(a) Even Functions (b) Odd Functions 

Figure 15. Index functions for the ID supercritical NLS equation 
computed at five values of 2.3 < a < 6.3, inclusive. 




(a) (b) 



Figure 16. The Kj as functions of cr for ID supercritical NLS. 

The Kj'^^ have the values indicated in Figure Moreover, 

kI""^ < 0, for 2.3 < cr < (74 
K^2^ < 0, for 2.3 < cr < 6.3 

where 

(4.29) 0-4 = 3.49928679909 

(e) (e) 

Proposition 4.11. Let Zl and , solve 
(4.30a) &^Z^'^ = lR + xB' 

(4.30b) &^Z['^ = 01 
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-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 
-0.7 



4 
(a) 




4 

(b) 




Figure 17. (4.31) as functions of for ID supercritical NLS. 



and define: 
(4.31a) 
(4.31b) 
(4.31c) 

The jj""* have the values indicated in 



£(e) ^(e) ^(e) 



7-(e) _ / /-(e) 7(e) 7(e) 



j: 



Figure | j7| 



Moreover, 



(4.32a) 
(4.32b) 

where 

(4.33a) 
(4.33b) 



4"'^ < 0, /or 0-5 < < 6.3, 
{4^ht^ - (ji'^V)/^^'^ < 0, for a,<a< 6.3, 



0-5 = 2.45649878, 
(76 2.45 3 79561. 



Proposition 4.12. Let u[°'' solve 
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(4.34) 

and define 

(4.35) K 

k[°^ has the values indicated in Figure 
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Moreover, 



K'f^ < 0, for 2.3 < cr < 0-7 
(77 = 6.1288520139 



where 
(4.36) 



4.3.3. Proof of the Spectral Property. To prove the spectral property in this case, we 
can use the orthogonahty of / to (t>2 to induce positivity of for aU 2.3 < a < 6.3, 
For orthogonahty to -R + xR' and (pi induces 



via Proposition 



4.10 



positivity for a-^ < cr < 6.3 through Proposition 4.11[ Finally, Proposition 4.12 
yields positivity of B'f ^ for 2.3 < a < a-j. 
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5. Discussion 

We have successfully proven there are no purely imaginary eigenvalues, both 
in the spectral gap and embedded in the essential spectrum, for a collection of a 
orbitally unstable solitary wave solutions of NLS. However, there is still much to be 
done. First, the range of orbitally unstable solitary waves has not been exhausted. 
We know there are no purely imaginary eigenvalues for cr > 2 in the ID NLS 
equation, and we anticipate that there are none for cr > 2/3 in 3D. We expect 
similar results for the entire unstable branch of CQNLS. Second, we have shown 



that the threshold found in |6], given by (1.10), is suboptimal; our range, (1.12|, is 
larger, but it too is likely suboptimal. 

This work also raises questions of what is satisfactory for a rigorous proof. We 
relied on the computer in the following steps: 

• Computing the solitary waves for the 3D, a nonlinear boundary value prob- 
lem on an unbounded domain; 

• Computing a finite number of index functions, linear initial value problems 
on unbounded domains; 

• Solving a finite collection of linear elliptic problems on unbounded domains; 

• Computing inner products on an unbounded intervals. 

To deal with these obstacles we truncated the intervals to [0, rmax] ([0, Smax] in ID), 
and introduced a finite number of discretization points within the intervals. Even 
with these computations, we have only established the results for isolated values 
of cr or 7; we argue by continuity that the results should extend to continuous 
intervals. 

However, we are able to produce a variety of consistency and a posteriori checks, 
such as verifying that the artificial boundary conditions are being satisfied suffi- 
ciently well. Despite our confidence in our computations, it would still be desirable 
to avoid the use of the computer. Indeed, an alternative method may permit us to 
move beyond the seemingly artificial restrictions on cr and 7 of our results. 
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Appendix A. Numerical Methods 

The codes used to produce the results appearing in this work are available at 
http : //www ■ math . toronto . edu/ simpson/f iles/ spec_prop_asad_simpson_cod^ . 

zip for examination, verification, and experimentation. 

The numerical methods we use are the same as those used in [12j , and separately 



in Simpson & Zwiers, 20 . For the 3D problems, we use the Fortran 90/95 based 
boundary value problem solver of Shampine, Muir, & Xu, fl9'. For the ID problems, 
we found it sufficient to use the closely related bvp4c routine in Matlab , [17^(18] , 
These algorithms require us to specify the problems as first order systems of the 
form 

where 5' is a constant matrix and f contains all non-singular terms, both linear and 
nonlinear. 

This automatically accommodates problems involving type singularities at 
the origin, as found in the equation for the 3D solitons and the k = harmonic 
boundary value problems. For higher harmonics, we use a point transformation 
W{r) — r''W{r) to transform 



Id k{k + d-2) 
dr 



W 



dr^ 



l + 2k d 
r dr 



W. 



A.l. Index Functions. For the 3D problems, we computed the index functions 
on the domain [0,200] with a tolerance setting of 10"^'^. The ID problems were 
computed on the domain [0, 100] with absolute and relative tolerances of 10"^'^. 
While Figures [3j [9| and [15] show the indicated numbers of roots, they obviously do 
not prove that there is no root at some larger value of x or r. 

To mitigate this uncertainly, there are two things we can check. First, we can 
examine the solution on a much larger domain. A more systematic approach is 
to note that since the potentials are highly localized, asymptotically, the index 
functions satisfy the estimates 

(A.la) f/C^^W « C^'V^ + C^'V^-"^-^ 

(A.lb) zW(r) « I^^^'V^ + i^J'^V^-^-^- 

where cj'^'' and D^'^'' are constants. Estimating these constants, and seeing that 
they have the correct signs and magnitudes, indicates that there cannot be further 
roots. For example, in Figure 20 we see that C^^'' ~ —0.3668 and C'f'' ~ —0.2393. 
The only root of 

-0.3668 - 0.2393/r 

is near —0.65. Thus, there should be no zeros in the "far field". Similar analysis 
works for the D^^^ constants. 



A. 2. Boundary Value Problems & Artificial Boundary Conditions. Since 
this algorithm requires us to compute on a finite domain, [0, r^ax]: wc must in- 
troduce artificial boundary conditions at r^^a^- These are discussed in the above 
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Figure 20. The index functions and the asymptotic constants for 
the 7 = .01 case of CQNLS in harmonic k = 0. 




oo o o o o 



"'*itf^i!iSV°"«.,a«- * » iS " 

40 



60 80 



100 



Figure 21. Verification of the artificial boundary conditions for 
the sohton, R, the eigenstate, 4>, and the solutions of the boundary 
value problems C/j"-* for the 7 = 0.01 case of CQNLS. 



references, but include 



9c^U=l^(rmax) + 9„L=1-R'(f'max) = 



for the 3D problems. Again, additional details are given in 12 20 . An example 
of an a posteriori check of these being satisfied is given in Figure |21[ which in- 
cludes computations of the soliton, the unstable eigenstate, and the solutions of 



the boundary value problems for Proposition 4.6 

The boundary value problems were solved on the domain [0, 100] in all cases. For 
the 3D problems, the tolerance was set at 10~^^. for the ID problems, the absolute 
tolerance was 10~^ and relative tolerance was 10^ 

A. 3. Root Finding. To compute roots of the inner products, or combinations of 
them, we wrap our routines within the fzero routine, with default settings, of 
Matlab ; we define a function that, for a given value of a or 7, returns the inner 
product, and pass this to fzero. We found this could provide approximately six 
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6 
4' 
2 



-2 



X 10 



-6^ 



4 6 

a - 0.8076980 



X 10 



Figure 22. Zooming in on the root of jf^ inner product of the 
3D NLS equation, the solution is appears piecewise constant. This 
is a closeup of Figure [6]Ja) . 



digits of accuracy in our 3D computations, notably for the bounds ( 1.12 ) and ( 1.13 1, 
while for the ID computations, we were able to get eleven digits. 

There appears to be a loss of sensitivity to the input cr or 7 for the 3D problems. 



near the root, j'f^ is piecewise constant. Refinements to our approach, with higher 
order corrections to the artificial boundary conditions, may improve the accuracy. 



An example of this appears in Figure 22 where we see that for a range of a values 
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